New Experimental Constraints on Non-Newtonian Forces below 100 fxm 
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We have searched for large deviations from Newtonian gravity by means of a finite-frequency 
microcantilever-based experiment. Our data eliminate from consideration mechanisms of deviation 
that posit strengths ~ 10'* times Newtonian gravity at length scales of 20 /im. This measurement is 
3 orders of magnitude more sensitive than others that provide constraints at similar length scales. 
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Physics beyond the standard model will be explored 
at the Large Hadron Collider in a few years, but gravi- 
tational measurements at small length scales are capable 
of investigating this region now. These measurements 
become particularly important in light of recent theoret- 
ical developments. Moduli, massive scalar particles that 
could mediate forces with strengths possibly ~ 10^ times 
Newtonian gravity, could have Compton wavelengths in 
the range of 10 fim-1 mm. These particles would mani- 
fest themselves as strong gravitylike forces. The hierar- 
chy problem, the seeming disparity (16 orders of magni- 
tude) between the standard model energy scale and the 
Planck scale, suggests other reasons to study this regime. 
This problem can be recast by asking why gravity is so 
weak compared to the other known forces in nature. Re- 
cent theoretical work P] suggests that extra spatial di- 
mensions, possibly as large as 1 mm and accessible only 
to gravitons, may supply a solution. Because of spread- 
ing into the extra dimensions, gravity would be diluted 
compared to the other forces, and therefore would seem 
quite weak. As gravity has not been well tested below 
the centimeter scale until recently 0, |^ , there are many 
reasons to investigate gravitylike forces at length scales 
of < 100 fim. 

Exotic massive scalar particles would modify the stan- 
dard Newtonian potential with an additional Yukawa- 
type term Q. A similar correction would exist for ob- 
servation of gravity at length scales roughly the size of 
any compactified extra dimensions Q. This leads to the 
following equation for the gravitational potential in the 
presence of such non-Newtonian effects: 
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Here G is Newton's constant, mi and m2 are the masses, 
r is the center of mass separation, and a is the strength 
(relative to Newtonian) of any new effect with a length 
scale of A. Non-Newtonian effects are typically parame- 
terized in terms of a and A. 

In this Letter we present data obtained for attracting- 
mass surface separations down to 25 /im, enabling inves- 
tigation of interaction scales below 10 /xm. 

The experiment is reminiscent of torsion pendulum 
gravity measurements, though an ac drive force is em- 
ployed, and the masses are much smaller. In place of 



the typical torsional fiber force sensor, a microcantilever 
was used. A mass attached to this cantilever was sub- 
jected to a time-varying force of a gravitational type, 
and the deflection of the cantilever (measured using fiber- 
optic interferometry 5]) provided a measure of this force. 
The experiment was performed at low temperature and 
in cryogenic vacuum to exploit the high quality factors 
attainable in microcantilevers under these conditions. 

The force sensors for these experiments were single- 
crystal silicon diving-board-shaped oscillators 250 ^m in 
length, 50 /im in width, and 0.335 /xm thick. The can- 
tilevers were fabricated using standard micromachining 
techniques and have quality factors in the 10"^ — lO'^ 
range (in vacuum). Spring constants of these sensors 
were 5.0-5.5 mN/m, and the resonant frequency of the 
lowest flexural mode of the cantilever was shifted from 
^ 7000 Hz to ~ 300 Hz upon mass attachment. 

Microcantilever sensitivity is usually limited by ther- 
mal noise, i.e., Brownian motion of the beam. This noise 
limit, analogous to the Johnson noise in a resistor, can 
be used to obtain the minimum detectable force of a can- 
tilever: Fmin — AkBTkB /ijJqQ. Here is the Boltz- 
mann constant, T is the temperature, k is the spring con- 
stant, LjQ is the resonant frequency (rad/s), Q is the qual- 
ity factor, and B is the measurement bandwidth. The 
typical force sensitivity for these (thermal noise limited) 
measurements was ~ 1 x 10~^^ N/\/Hz. With averaging 
times of a few hours, the present devices could reach an 
ultimate sensitivity below '--^ 1 x 10~^^ N (200 times the 
expected Newtonian gravitational force) . 

The force measured was that between a mass placed on 
the end of the cantilever (the test mass) and a larger mass 
oscillated a small distance away (the drive mass). The 
test masses were gold rectangular prisms, 50 x 50 x 30 fiu^ 
in size and ~ 1.4 /xg in mass, and were attached to the 
cantilevers with a thin (typically ^ 1 /im) layer of epoxy. 
The drive mass was constructed of ten bars: five gold 
bars alternating with five silicon bars. Each bar was 
100 /im X 100 /tm x 1 mm long. (See Fig.^for schematic 
representation.) This construction, when oscillated in 
the plane of the drive mass and perpendicular to the di- 
rection of the long dimension of the bars, is expected 
to gravitationally excite the test mass at a harmonic of 
the oscillation frequency. This shift of the gravitational 
signal frequency with respect to that of the drive mass 
oscillation helped prevent spurious excitation of the can- 
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FIG. 1: Schematic top and side views of test mass on can- 
tilever and drive mass below, showing alternating gold and 
silicon drive mass bars. Drive mass motion was in the Y di- 
rection. 



tilever due to unwanted vibration at the signal frequency. 

To test the system and precisely align the drive mass 
with the test mass on the cantilever, as well as to pro- 
vide an in situ equivalent measurement, a magnetic ana- 
log of the gravitational experiment was used. The gold 
bars in the drive mass were connected electrically at al- 
ternating ends (see Fig. ^ to form a meander. When 
a do electric current is driven through the path defined 
by the gold bars and their interconnects, a space- varying 
magnetic field (with a spatial periodicity half that of the 
gravitational response) is created above the drive mass. 
This field couples to the test mass through an evaporated 
layer of nickel on the test-mass surface closest to the drive 
mass. The oscillation of the drive mass creates a time- 
varying magnetic field at the location of this layer. A 
signal proportional to current through the meander ver- 
ifies the system's functionality. 

To fabricate drive and test masses, bulk silicon was 
etched to form a mold into which gold was evaporated. 
After evaporation, the mold was ground and polished. In 
the case of the test masses, the silicon was then dissolved 
to release the gold prisms (for fabrication details, see Q)- 

The oscillation of the drive mass beneath the test mass 
was facilitated by attachment of the drive mass to the free 
end of a clamped piezoelectric bimorph actuator. The 
actuator was 1 mm thick and 4 mm wide and had a free 
length of 41 mm. 

The actuator's motion was characterized ex situ using a 
laser beam-bounce method. In the experimental probe, 
the actuator's motion was measured capacitively using 
a calibration obtained from this initial characterization. 
The actuator's drive frequency was tuned to one-third 
of the cantilever resonant frequency by means of mag- 
netic excitation. The amplitude of actuator motion was 
98±7 fim. Finite element calculations were performed to 
determine the expected Newtonian and non-Newtonian 
response of the test mass. 

The nonlinearity in the piezoelectric actuator caused 
a small amount (4%) of actuator vibration at the can- 
tilever resonance (the third harmonic of the actuator fre- 
quency) and could lead to a spurious signal. Therefore, 
two mass-spring vibration isolation stages separated the 
cantilever mount from the actuator mount. Each stage 



had a resonant frequency of '--^ 2 Hz for both vertical and 
horizontal motion, and thus together the stages provided 
an attenuation of ^ 10^ at 300 Hz between the mounts. 

A stiff metallized shield was placed between the drive 
and test masses to minimize the effects of electrostatic 
and Casimir forces by preventing direct ac coupling be- 
tween the masses. The shield was a 2.8 mm by 5.2 mm 
silicon-nitride membrane 3 /im thick onto both sides of 
which was evaporated 100 nm of gold. Supported by a 
1 cm^ silicon wafer die, the shield was attached to the 
cantilever wafer 15 /im from the cantilever. 

The vertical (Z) and in-plane {X, Y) relative displace- 
ments of the test and drive masses were determined by 
means of a capacitive position sensor (CPS) similar to 
that described in 0. Two simpler capacitors, in con- 
junction with the CPS, provided measurement of rela- 
tive tilt. After cooling, the position was adjusted until 
the original capacitive readings were regained, using a 
three-axis translation stage. Because of the drive mass 
geometry, the position in X (the direction along the drive 
mass bars) need only be determined to an accuracy of 
~ 100 /im. In the Y direction, greater precision was re- 
quired as the gravitational response was expected to vary 
by ~ 100% over 50 /im in Y. The CPS enabled deter- 
mination of relative position in this direction to ±2 /im. 
The tilt was adjusted by differential heating of the upper 
vibration isolation stage support springs. 

After lateral realignment, an absolute value of the ver- 
tical separation between the masses was determined by 
means of direct mechanical contact (sensed by the can- 
tilever) between the drive mass and the shield. Subse- 
quent motion away from the shield was measured with 
the CPS, and the separation between the masses was 
thus known to an accuracy of 2 /im. 

The cantilever, actuator, and vibration isolation stages 
were suspended in a vacuum can at the end of a probe 
that was inserted into a liquid helium research dewar. An 
exchange gas space separated the inner vacuum can from 
the liquid helium in order to reduce external vibrations 
due to helium boiloff. The entire system was supported 
by 2 m long ~ 2 Hz springs for additional isolation. The 
experiments were performed at T = 9-11 K, and at gas 
pressures of less than 10~^ torr (room temperature read- 
ing) with further cryopumping while cold. 

The sinusoidal actuator voltage and the cantilever dis- 
placement signal were sampled and stored simultaneously 
at 10 kHz using an analog-to-digital converter. The dis- 
placement data stream was then sorted into bins using 
the phase of the actuator drive voltage and averaged by 
actuator phase. As the actuator phase is a measure of the 
lateral position of the drive mass, this technique is similar 
to a lock-in technique in that signals that are not phase 
coherent with the actuator drive signal are rejected. 

The third harmonic of the Fourier transform of the av- 
eraged data is the motion of the cantilever that is phase 
correlated with the third harmonic of the drive mass mo- 
tion, i.e., the gravitational- like response of the test mass 
[see Fig. |2Ia)]. Using the spring constant and the mea- 
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FIG. 2; (a) Actuator voltage (left axis), raw cantilever dis- 
placement (right axis), and averaged cantilever displacement 
(right axis, 10 min averaging time) as a function of time over 
one period of actuator motion, (b) Averaged data for one pe- 
riod of actuator motion for five values of dc current through 
the drive mass meander. 



sured quality factor of the cantilever, the force on the 
test mass due to the drive mass can be determined. 

Data were collected as a function of drive mass me- 
ander current and Y position in order to align the drive 
mass to the test mass. Figure |2Ib) shows data for one 
period of the drive mass actuator motion as a function of 
drive mass meander current. The phase of the magnetic 
response changes by tt as the current is changed from 
negative to positive. Also present is a small asymmetry 
in the magnitude of the magnetic signal for positive and 
negative values of the nominally same current magnitude. 
The observed zero-current signal adds phasewise to the 
signal for a particular current, producing this asymmetry. 

The third harmonic of two sets of data acquired over 
1000 sec and averaged as a function of measurement time 
is plotted in Fig.|31(Q = 1180 for these data). Also plot- 
ted is a line representing the theoretical thermal noise of 
the cantilever. Data taken while the actuator was oscil- 
lating far from the test mass match the theoretical curve 
within error and imply that the dominant experimen- 
tal noise is thermal in nature. (Note that the theoretical 
curve depends on experimental parameters.) Data taken 
with the drive mass oscillating and in close proximity to 
the test mass show a signal that clearly is not dimin- 
ished with averaging and is patently above the thermal 
noise. According to finite-element calculations yielding 
an expected gravitational response of 4.0x10"^° N, this 
anomalous signal is too large (by a factor of ~ 2000) to 
be of Newtonian origin. Within experimental error, the 
magnitude of the signal (relative to the actuator drive) 
as a function of F-direction offset is not consistent with 
a mass-dependent force. 

These data lead to the conclusion that the experi- 
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FIG. 3: Background and experimental signal, as well as the- 
oretical thermal noise, as a function of averaging time. The 
experimental signal levels off at a force of 8.9 x 10"^'^ N. 



ment's force resolution is presently limited by an envi- 
ronmental effect, most likely an electrostatic interaction 
between the shield and the cantilever as discussed be- 
low. While the shield prevents direct coupling of this 
sort between the masses, a secondary effect due to shield 
motion could spuriously drive the cantilever. Other pos- 
sible non-mass-dependent effects are too small to appear 
at this level. For instance, randomly quenched magnetic 
impurities in the gold of the drive mass produce a mo- 
ment that would result in a force below the level of New- 
tonian gravity, and coordinated pressure variations in the 
ballistic, residual gas in the probe are similarly small. 

Shield motion was characterized (using the fiber inter- 
ferometer) under standard experimental conditions. The 
response was not sensitive to the current through the 
drive mass meander or relative mass position. With the 
actuator swinging its full amplitude, the shield motion at 
the first harmonic was ^ 10 pm. The signal at the third 
harmonic was an order of magnitude less. 

Motion of the membrane could drive the cantilever by 
creating an ac Casimir force. An upper bound can be 
calculated assuming the cantilever is a conductor and 
using the Casimir result |lC| for two conducting planes 
(force per unit area), Fa = — tt"^ he/ 240 z'^. For 10 pm of 
motion of the shield, the oscillatory force is on the order 
of 1 X lO"^'^ N, much less than the thermal noise of the 
cantilever for experimental measurement times. 

The motion of the membrane could also induce can- 
tilever motion if there were a significant potential dif- 
ference between the metallized top of the shield and the 
cantilever. Making a worst-case assumption of two paral- 
lel plates with an area the size of the cantilever, a voltage 
difference 0.3 V would produce a force on the cantilever 
of 1 X 10"^^ N. Potentials of this size have been observed 
in other experiments jlll | , and we cannot yet exclude the 
possibility that such a force currently limits our experi- 
mental sensitivity. Shield motion was most likely due to 



4 



r-TTi — 1 

r ^\ 
I gaugs\ 

r bosonX 

- -i \ \^ 

r \ 

r \ ^ ^ 

r 


— 1 1 1 1 1 1 1 1 1 

Lamoreaux 1 

Region excluded 
V^^^ by experiment ■ 

V UColorado 

\\* i 


1 rmocluli 


J\~ \- ; 


|. dilaton 

r 


1 \^^^/ U W ashington 1 


r This work 




r 


Irvine 1 


r 





A,(|um) 

FIG. 4: Strength versus length scale parameter space [see 
Eq. Q] for non-Newtonian effects showing area excluded by 
present (darker shaded region bounded by signal plus 2a- 
level error) and previous (lighter shaded region) experiments. 
Lines labeled Lamoreaux, UWashington, UColorado, and 
Irvine are from [T^, Q, E3i ^^'^ H^ . respectively. Theoreti- 
cal predictions (dashed lines) are adapted from |4l ll5|l (dilaton 
and moduli) and from (gauge bosons). 

inhomogeneities in the drive mass surface caused by the 
use of electrically dissimilar materials. The design of up- 
coming experiments contains insulating and conducting 



layers above the drive mass bars; the conducting layer 
will shield the effect of any drive mass surface variations. 

The existence of a spurious force of average magnitude 
8.9 X 10""'^^ N allows calculation of an upper limit on the 
magnitude of any gravitylike force. Figure 01 displays the 
region of non-Newtonian interaction strength versus in- 
teraction distance parameter space that can be excluded 
by these results jl7l |. as well as other current experimen- 
tal bounds. For a length scale of 15 /im, the present data 
improve constraints on the strength of non-Newtonian ef- 
fects by '--^ 3 orders of magnitude. The region eliminated 
is an important one, as it contains previously unexplored 
space for scalar moduli particles. 

We have presented results from a search for gravity- 
like deviations from Newtonian theory at distance scales 
below 100 /im using a microcantilever approach, with 
masses of size on the order of these scales. Our data of- 
fer a new limit on non-Newtonian effects in the range of 
10 /im, constraining theoretical predictions of moduli. 
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Note added. — After submission of this manuscript, re- 
sults for larger length scales were published from a similar 
experiment at the University of Colorado jl^] . 
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